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Abstract 

The motion of a rigid body immersed in an incompressible perfect fluid which occupies a three- 
dimensional bounded domain have been recently studied under its PDE formulation. In particular 
classical solutions have been shown to exist locally in time. In this note, following the celebrated 
result of Arnold [I] concerning the case of a perfect incompressible fluid alone, we prove that these 
classical solutions are the geodesies of a Riemannian manifold of infinite dimension, in the sense that 
they are the critical points of an action, which is the integral over time of the total kinetic energy of 
the fluid-rigid body system. 
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1 Introduction 

We consider the motion of a rigid body immersed in an incompressible homogeneous perfect fluid, so that 
the system fluid-rigid body occupies a smooth open and bounded domain C M'^. The solid is supposed 
to occupy at each instant t ^ a, smooth closed connected subset S{t) C which is surrounded by a 
perfect incompressible fluid filling the domain T{t) := il \ S{t). 

For the point of view of PDEs, this system have been recently studied in [S], [S], [TU], [5], [S] which 
have set a Cauchy theory for classical solutions. 

The aim of this note is to provide a rigorous proof that the classical solutions can be equivalently 
thought as geodesies of a Riemannian manifold of infinite dimension, in the sense that they are the 
critical points of an action, which is the integral over time of the total kinetic energy of the fluid-rigid 
body system. It was pointed out in a famous paper by Arnold [1] that both the Euler equations for 
a rigid body as well as the Euler equations for a perfect fluid can be derived with this approach. The 
motion of a rigid body in a frame attached to its center of mass can be considered as a geodesic on the 
special orthogonal group S0{3). On the other hand the motion of a perfect fluid filling a container fl 
(without any immersed rigid body in it) can be considered as a geodesic equation on the space Sdifr'"(57) 
of the volume and orientation preserving diffeormorphisms of ft. 

It is hence natural to try to extend this analysis to a system of interaction of a perfect fluid and a 
rigid body. In particular cases, when the fluid is irrotational or when the vorticity of the fluid is given by 
a finite number of point vortices, so that the dynamics is finite-dimensional, this was studied in details in 
pn [T^ ; see also references therein. The goal of this paper is to observe that one can see the motion of a 
rigid body in a fluid governed by the incompressible Euler equations, as a geodesic flow, in the presence 
of a regular distributed vorticity, as well. 
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The structure of this paper is as follows. In Subsection 11.11 we first recall the PDE formulation of 
the system. Then in Subsection ll.2l we describe the infinite-dimensional manifold and the action used in 
the geometric formulation of the problem. In Subsection 1 1.31 we state the main result of the paper, that 
is, the equivalence of the two points of view. Section [2] contains the proofs of the various statements. 
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1.1 PDE formulation 

The dynamics of this system can be described thanks to the following PDEs: 

(1) 
(2) 
(3) 

(4) 

(5) 
(6) 

(7) 
(8) 
(9) 
(10) 

Equations (H]) and ^ are the incompressible Euler equations. The vector field u is the fluid velocity and 
the scalar field p denotes the pressure. Equations ([3]) and (|1]) are Newton's laws for linear and angular 
momenta of the body under the influence of the pressure force. Here we denote by m the mass of the 
rigid body (normalized in order that the density of the fluid is pp = 1), by xsit) the position of its 
center of mass, n{t, x) denotes the unit normal vector pointing outside the fluid and dT(t) denotes the 
surface measure on dS(t). The time-dependent vector 

m x'sit), (11) 

denotes the velocity of the center of mass of the solid and r denotes its angular speed, so that the solid 
velocity field is given by 

v{t, x) e{t) + r{t) x{x- XB{t)). (12) 

In (U]) the 3x3 matrix J = J{t) denotes the moment of inertia which depends on time according to 
Sylvester's law 

J = QJoQ*, (13) 

where Jq is the initial value of J" and where the rotation matrix Q e 50(3) is deduced from r by the 
following differential equation (where we use the convention to consider the operator r(t) x • as a matrix): 

Q'{t) = r{t) X Q{t) and Q{0) = Ida • (14) 

The matrix Jb can be obtained as follows. Given a positive function pg^ e L°°{So;R) describing the 
density in the solid (again normalized in order that the density of the fluid is pp — 1), the data m, xq 
and Jo can be computed by it first moments 

m :— / psgdx > 0, (15) 
J So 

mxo := / xpsg{x)dx, (16) 
Jso 
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Jo / pso{x){\x - xo\'^Id3~{x - xo) {x - xo))dx. (17) 

Finally, the domains occupied by the sohd and the fluid are given by 

Sit) = ^^^XB{t) + Q{t){x-xo), X e 5o| and =n\'S(J), (18) 

starting from a given initial position Sq C such that J-q :~ fl\So. Let us underline that S{t) and dQ 
being compact, and since S{t) C fl, the solutions that we consider satisfy 

dist{S{t),dn) > on [0, T]. (19) 

Let us give a precise definition of the classical solutions examined in this paper. 

Definition 1 (Classical solutions). We call classical solution of the PDEs formulation on [0,T] some 

(for some A G (0, 1)) satisfying ((H -((H]) . 

The local-in-time existence and uniqueness of classical solutions to the problem ([T|)- P^ holds when 
the initial velocity of the fluid is in the Holder space C^''' cf. [5] for a precise statement. Let us also 
mention the earlier results of Ortega, Rosier and Takahashi where the body- fluid system occupies 

the plane M^, Rosier and Rosier [TU] in the case of a body in R'^ and Houot, San Martin and Tucsnak 
[B] in the case (considered here) of a bounded domain, with the initial velocity in a Sobolev space H™, 
m ^ 3. 

1.2 Geodesic formulation 

Let us now turn to the geometric viewpoint. We flrst describe below the inflnite-dimensional space of 
configuration of the system. Next we introduce a natural action, which allows to define our notion of 
geodesic. 

1.2.1 Rigid movements 

Let us first describe the rigid part of the motion. To a velocity vector field w of a rigid body one associates 
the flow 

dtT{t, x) = v{t, T{t, x)) and t(0, x) ^ x for {t, x) G [0, T] x Sq. (20) 
It is easy to integrate to find 

T{t,x) = Xsit) + Q{t){x - Xo), 

where 

XB{t) ■■= T{t,xa), 

and where Q(t) is obtained from v by (ITil) and r is given by p^ . 

The flow T can be seen as a function of the time with values in the Lie group SE{3) of rigid 
motions (the special Euclidean group), that is the group generated by translations and rotations in 3D. 
Its tangent space at Id G SE{3) is 

se(3) := ridS'£;(3) = G C1(R=';R3) / D{v) = o|, 

where D{v) denotes the tensor of deformations 2D{v) := [diVj -I- djVi)i,j. Given xb in M'^, we have the 
following: 

se(3) = |u : ^ / 3(£, r) G M'"^ x M? , Vx G M^ v{x) ^ £ + r x {x - xb)}- 

Moreover given xb in M.'^ and v G se(3), the ordered pair {£,r) above is unique. Hence se(3) x 
Skew3(R). 
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Accordingly, the tangent space of SE{3) at r G SE{3) is 

TrSE{3) := |w o T with V £ se(3)|. 

Now, xq being given, we introduce the following projections on TrSE{3): to a = v o t € TrSE{3) we 
associate {L-rla], Rr[a]) the unique ordered pair {£,r) associated to v with xb '■— t{xo), in other words: 

a{T^^{x)) — v{x) — Lr[(T] + i?r [c] X (a; — T{xoj). 

Let us conclude this subsection by describing the energy of the solid. Using the choice of xb (t) as 
the center of mass of the body at time t, we have that 

PSoirt-\x)){x~XBmdx^O, (21) 

Sit) 

and therefore, for any (^i, ri, i!2, ^'2) G (M^)'*, for any t, 

pSo{Tt'\x)) {£1 + n X (x - XB{t))) ■ {£2 + r2 X (a; - x^W)) dx = m£i • £2 + J{t)ri ■ rz, (22) 

s{t) 

where J{t) is given by and the notation r^"^ stands for the inverse of the function Tt :— T{t, •). 

1.2.2 Fluid displacements and Arnold's geodesic interpretation 

Let us briefly recall Arnold's interpretation of the Euler equation. To a velocity vector field u satisfying 
the incompressible Euler equations in ft (without body) one associates the flow rj defined on [0,T] x £7 

by 

dtrjit, x) — u{t, ri{t, x)) and 7^(0, x) — x. (23) 

The flow 7] can be seen as a continuous function of the time with values in the space SdifF''(fi) of the 
volume and orientation preserving diffeormorphisms defined of J7. The latter is viewed as an infinite- 
dimensional manifold with the metric inherited from the embedding in L^(0.] 1^^), and the tangent space 
in ry e SdifF^(r2) is 

T,, Sdifl:+(r2) := |mo7; with ii e C^(rj;R^) such that div (u) = in Vt and u-n = on arjj. 

Euler's equations are then interpreted as a geodesic equation on Sdifr'"(51). The pressure field appears 
as a Lagrange multiplier for the divergence-free constraint on the velocity. Ebin and Marsden proved in 
[3] the existence of these geodesies when the initial velocity of the fluid is in the Sobolev space H^{VL), 

1.2.3 Possible configurations of the fluid/body system as a Riemannian manifold 

To introduce the geodesic formulation of the motion of a rigid body immersed in an incompressible perfect 
fluid, we flrst describe the infinite-dimensional manifold on which these geodesies will be considered. 

The set of the possible configurations as a Riemannian manifold. To begin with, we first 
describe the set of the possible configurations of the system at a fixed time by setting 

C := I (r,77) e SE{?,) x C^^^{Fo;M?) such that t(5o) C n, 

77 is a volume and orientation preserving diffeomorphism J-q — >■ i7 \ [r(5o)] |. 

We will represent (r, 77) by : £7 — 51 such that ~ r and (t)\jr^ — rj. Note that (p is not necessarily 
continuous. 



4 



Let us observe that according to the Helmoltz decomposition, in C^'^(Jo; M^), the space of divergence- 
free vector fields is closed and admits a topological complement. Therefore, one can see that C is a 
submanifold of a manifold modelled on the Banach space 

E :=se(3) x C^'^{TQ;m.^). 

We will be interested in the tangent space of this infinite-dimensional manifold. Let us first recall that by 
definition the tangent space T^^-^jC in (r, 77) G C is the set of equivalence classes of germs of continuously 
differentiable function 9 : (— e,e) — !• C (for some e > 0) such that 0(0) = (t, 77), for the equivalence 
relation: 

61 - 62 iff e'i(O) = 62(0) in local charts. 

Equivalently (cr, /i) £ T^^ ^j-^C is the set of equivalence classes of pairs {ip, e) where "0 is a chart defined on 
a neighborhood of (r, r/) E C and e G £^ for the relation 

(V'i,ei) = (^'2,62) if D{ip2 0iJ;^^){ipi{T,r])) ■ ei = 62. 
Clearly this makes T^^jj^C a linear space and TC := (J {{{'''^'n)} ^ -^(r,?;)^^ a vector bundle whose 
base space is C. 

Let us introduce the following notation: given (t, 77) G C and fi G C^''''(J^o; we introduce Ujj[ii] : 
rji^o) ^ M3 by 

Recall that due to the openness of ft and closedness of Sq one has d{(j){So), 9^) > for G C. 

Now the tangent space of C is described by the following proposition. 
Proposition 1. Let {t,vi) G C. Then using the notations 

XB := t(xo), S := r(5o) and T:^n\S, (24) 

we have 

T>3\ I 



T(,,^)C = I (a. A*) G T,SE{i) x C^-'^{T^- 



div([/^[^]) =0 m JT, (25) 
Ur\[i\{x) ■ n{x) = {Lr\G\^ Rr\a\'%. {x — xb)) ■n{x) on dS, (26) 

L/^[^](a;) ■n(x) =0 on dil^. (27) 

before, n denotes the normal unit vector on dfl and dS , pointing outside the fluid. 

The proof of Proposition[T]will be given in the next section. We will represent {a, fi) hy 14 o cj) : —i' M.^ 
where : 57 ^> M'^ is given hy U^gor — a and U\jr = ?7,j[/i]. In the same way as 0, U can be discontinuous. 

The manifold C can be endowed by the following Riemannian metric: for any (j) E C, for any Ui o 0, 

U20 (l> & T^C, 

Iplio (j),U2 (j))^ := I po{(p^^)Ui ■U2dx. 
Jn 

Above we used the density po defined on by: 

Po{x) -- 

Splitting the fluid flow and the solid one this reads: for any (ai, /ii), (cr2, M2) in T'(t.?))C, we have 



1 in J^o, 
PSoix) in Sq. 



{{ai,fii),{(72,tJ-2))(T.r,)^mLr[a-i]-Lr[a-2]+J[T]Rr[(Ti]-RT[o-2]+ / [/^i] • CX,, [^2] da;, (28) 
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where J\t\ is the inertia matrix deduced from the initial one by the rigid transformation r that is 

J{A = Q[r]JoQ[r]*, (29) 
where Q[t] is the rotation matrix eanonically associated to r, that is, its hnear part. 

Let us stress that this metric defines a weaker topology than the original one on C. 

Curves of configurations. We now turn to time-dependent displacements. Let us be given (to,??o) 
and (ri,?7i) in C. We introduce 

£ := I (r,r/) £ C\[0,T];C), such that: 
i. t(0) = To, ?7(0) = T]o, 

it. r{T)^Ti, v{T)^m }■ (30) 

It is easy to verify that £ is a submanifold of a manifold modelled on the Banach SpcLC6 Ej" '. — 
CH[0,T];se(3)) x C^'^{[0,T] x J"o;K^)- 

The tangent space of this manifold is described by the following proposition. 

Proposition 2. Let {r^rj) e C. Then using the notations 

XB{t) := Tt{xo), S{t) := Tt(5o) and F{t) := n \ S{t), (31) 

we have 

T^r^.^)C = {(a,/i) e CWQ,T]-Tr,SE{i)) x C^'^{[Q,T] x F^-.W") j 

cr(0) = f7(T) = and fi{0) = ^(T) = 0, (32) 
div {U^ifi]) = in Fit) for each t G [0, T], (33) 
U^[^i\ ■ n = {Lr[a-] + Rr[cT] x [a; - XB{t)]) ■ n on dS{t) for each t £ [0,T], (34) 

L/J/i] ■ n = ondfl for each i G [0, T]|. (35) 

Remark 1. Here we make the abuse of notations C^([0,T];Tr^SE{3)), since Tr^SE{3) actually depends 
on t. We consider a as the section of a vector bundle rather than as a function. One can for instance in- 
terpret {[0,T];Tr^S E {3)) as the set of mappings a such that {t,x) i-^ a{t,Tf^{x)) is in C^{[0,T];sc{3)). 
Also, we sometimes drop the dependence of the objects on t in order to simplify the notations. 

Proposition [5] is the time-dependent counterpart of Proposition [TJ We do not provide a proof since it is 
only a matter of adapting the proof of Proposition [T] with a harmless parameter. The only new point is 
to observe that the extremities of the curves being prescribed (the conditions i) and ii) in (pOU the fields 
(7 and fi vanish when i = or T. 

As previously we will represent (r, 77) G £ by (/> : [0, T] x J7 ^ 57 such that for any t G [0, T], 

Ht, ■)\so = Tit, •) and <j){t, Obo = v{t, •)• 

We represent (cr,^) G T(^r.r,)'C- by W o : [0,r] x ^2 ^ ^^gre U : [0,T] x 17 gyyen, for any 

tG[0,T],by 

U{t,-)\sit)OT{t,-)='j{t,-) and U{t,-)\T{t)^U^Mt,-). 

1.2.4 The geodesic interpretation of the motion of a rigid body immersed in an incom- 
pressible perfect fluid 

Here we consider geodesies as critical points of the following action on the manifold C: 



Ai(j)) ■.= \ I po\dtc^\^dxdt. (36) 
^ -'[o,T]xn 
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We see that the action is obtained by integrating the squared norm (associated to the metric of C) of 
the tangent vector to the curve 0, that is 

-4(0) = 11 {dtcbit,-),dt(bit,-)}4,{t,-)dt. 

^ J[0,T] 

Separating the fluid and the body parts in the integral, and using (P^ . we see that 

A(t>) = ^ / (mm\^ + Jl-ritMt) ■ rit) + f \u{t, x)\'' dx) dt, 

where 

i Lr{dtT\, r -.^ Rr{dtT\ and w [9^77]. (37) 

In this writing, we recognize the integral over time of the kinetic energy of the fluid-body system. 

Moreover, going back to p6p . since the action is a continuous quadratic form on £, we deduce that 
A is differentiable on L with 

Z?^(0) • (iY o 0) = I podt(t>-dt{Uo(l))dxdt= [ {dt(t>{t,-),dt{U o <t>){t,-))^(^t.-)dt. (38) 

"'[O.TJxn "'[0,T] 

This leads us to the following natural definition. 

Definition 2. We say that cf) G C is a geodesic on C if for any U o (p £ T^C, DA{(j)) ■ {U o (p) — 0. 

1.3 Equivalence of the two points of view 

The main result of this paper is the following. 

Theorem 1. If {u,XB,r) is a classical solution of the PDEs formulation on [0,r] then {T,ri) defined by 
formulas (|20p and (j23p is a geodesic on C. 

Conversely, let {t,7]) C he a geodesic. Then {u,XB,r) where {l,r,u) is defined by (|37|) and xb 
is obtained by pi?]) and pT|). is a classical solution of the PDEs formulation on [0,T] (in the sense of 
Definitions^. 

Theorem [T] will be proved in Subsection 12.21 



Remark 2. Let us define the length of a curve 4> € C: 

A((/.) := / {{dt<P{t,-),M{t,-))^^t^.))Ut, 

J[0,T] 

and consider 

d := inf A((/)) 

^inf/ l\iit)\^ +J{t)rit)-r{t)+ [ \u{t, x)\'' dx] dt, (39) 
Jlo,T] \ J-ntiTo) ) 

where the infimum is performed over (j) = (r, r/) G £. We should say that d is the geodesic distance between 
the configurations {Ta,rj^) and [Ti,r]i) ofC. If {T,ri) S C realizes this infimum and is parametrized by t 
in such a way that the energy does not depend on time then (r, 77) also minimizes the action A over C. 
Conversely, by the conservation of energy, any geodesic is parameterized proportionnaly to arc length. 

Let us mention here two open problems. 

Open Problem 1. Is it possible to prove that for T small enough and (t, 77) G C such that the 
associated {u, xb, r) is a classical solution of the PDEs formulation on [0, T], one has for any (f, fj) G C, 
A{t, 77) ^ A(t, fj), with equality if and only if (f , 77) = (r, 77)? This should extend the result obtained by 
Brenier cf. [2] in the case of a fluid without body. 
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Open Problem 2. Is it possible to adapt the strategy that Ebin and Marsden used in [3] in the case 
of a fluid alone to the case with a body, that is to prove the existence of a free torsion connexion and 
of some geodesies by using parallel transport, and then to prove that these solutions also solve the PDE 
formulation? 

Let us also mention studies connected with the stability properties of the system. In the case of a fluid 
alone, Arnold T uses a notion of Riemannian curvature to investigate the stability of two-dimensional 
stationary flows. In the case considered here of a fluid-body system, this stability was studied by Ilin 
and Vladimirov [131 IE] . 



2 Proofs 

In this section, we prove the claims of Section [TJ The existence of a solution of the PDE system will not 
be needed. 

2.1 Proof of Proposition [T] 

1. Let us first show that any {a,fi) satisfying (|25p - (P7)) belongs to Tfj-^^-^C. First, we define 9s £ 
C\{-e,e)-SE{i))hy 

dsOsisM = vs{s,es{s)[x\), esmx] = t[xI 

where 

vs{s,x) := Lr[cr] + Rt[(j] x (x - ^^(s)), xb(s) xb + sLr[a], 
with £ > small; in particular we can require that for all s £ (— one has 

Ss ■.= esis)[s] c n. 

Lemma 1. There exists a family ^pg : J- ^ J-g :^ fl\ Sg, for s £ (— e, s), smooth in its arguments, such 
that 

V-o = Id^, (40) 

Vs G (— e,£), = ^5(5) ^'^ neighborhood ofS, (41) 

Vs £ (—£,£), ips = Id in a neighborhood of dfl, (42) 

Vs e (-£,£), Vx e J", det(V:rV's) = 1- (43) 

Proof. We define 

V{s,x) := -^cuT\{cj,{x){x X + \\x - XB{s)f Rr[a])) , (44) 

where is a smooth function equal to 1 in a neighborhood of S and in a neighborhood of dfl. We 
associate ips as the solution of 

dsil^s{x) = V{s,'ips{x)), ipoix) = X. 



Equations (|40l) and (|42)) are straightforward. Equation (|4T|) follows from the fact that V coincides with 
Vs in some neighborhood of dS, reducing e if necessary. Finally, since V is clearly divergence-free, P5)) 
follows from Liouville's theorem. □ 

Now we construct the vector field: 

vr{s, ■) := F(s, •) + {^s)*{UM - ^(0- •)) in J".- (45) 

(Recall that the push-forward is defined as {ip^,v){y) = d^^-i i^y-f[v{^~^ (y))].) 

Now we introduce a corresponding fiow 6jr according to the s-variable, starting from rj: 

ds9jr[s, X) = Vjr[s, 9jr{s, x)) , lu (-£, £) X J'q , 

6jr{0,x) = r]{x) in Fq. 
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(One can for instance smoothly extend w jr(s, •) to R" to define this flow in a standard way.) 

It remains to check that 6 {0s,6^) e C^{{-e,e);C), that 9(0) = (r, 77) and that e'(0) = {ct,^i). 
The two latter claims come directly from the construction. 

Let us prove that 9jr is volume-preserving. To that purpose, we first notice that div (U^) — div (V{s)) = 
0. Using the fact that the push- forward of a divergence- free vector field by a diffeomorphism with unit Ja- 
cobian determinant is still divergence-free (see for instance [71 Proposition 2.4]), P5)) and pS)) . we deduce 
that vjr(s) is divergence-free. Hence it follows that 9jr is volume-preserving by Liouville's theorem. 

The main point, that is that, for any s G (— e,e), 9j^{s) sends to J> can be seen as follows. It 
suffices to prove that 

vjr{s,-).n ^ vs{s,-).n ondSs, (46) 
vjr{s,-).n = ondn. (47) 

Using (H^ and (HH), we see that vjr = Urj[p] in some neighborhood of dil, so (|T7)) is a consequence of 
(P7)l . For what concerns from we see that C/,,[/i] — ^^(O,-) is tangent to dS. It follows that 

{ips)*{U^[fi] — vs{0, •)) is tangent to i9<Ss, which gives 

The other requirements, namely that djr is orientation preserving and has the claimed regularity, are 
clearly satisfied. 

2. Reciprocally, given {<J,ii) in T^^j^jC, by definition of the tangent space, there is = {QstOj^) G 
C^{{~e,e)]C) such that 9(0) = {r.-q) and that 9'(0) = (ct,^). Let us check that are satisfied. 

First, it is obvious that a = 9g{0) G TrSE{3). Next one can define 

vjr{s,x) := ds9jr{s,9~^^{s,x)), 

as a vector field on J^g. Then as 

a,0^(s,x) =i;^(s,0^(s,x)), (48) 

and since 9jr is volume-preserving, by Liouville's theorem we infer that div vjr = 0. Using again (j48p 
and the fact that 9jr sends Tq to Tg, we also easily infer (^5)) and ([Tf)). □ 



2.2 Proof of Theorem [T] 

We start with the following lemma. 

Lemma 2. For any (p ~ (r, rj) G C. for any U o <f) £ T^C, 
DA{(I)) • (W o 0) = 

/ (me{t) ■ dtLr[<y][t) + J[T]r{t) ■ dtRr[<j]{t) + [ dtr]{t,x) ■ dt{U{t,r]{t,x)))dx) dt 

where t and r are given by (|37p . Q[t] is the linear part of t and J^[t] is given by (j29p . 
Proof. Let us split the fluid flow and the solid one in ((55)) and prove that 

/ pso dtT ■ dta dx dt ^ / (me ■ dtLr[(j] + J[T]r ■ dtRric^]) . 

J[0,T]x5o JlO.T] ^ ' 

Now in order to prove (j49p we first perform a change of variable to get 

/ psodtT ■ dtadxdt^ / / psa{Ti'^{x))dtT{t,Tir'^{x))-dtCF{t,Tf^^{x))dxdt. (50) 

^[0,T]xSo J[0,T]Js[t) 

Now, by definition of e and r we have 

^tT{t,T^\x))^e + rx{x-XB)■ (51) 



(49) 
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In particular, taking the linear part of these afRne transformations, and using that the linear part of dtT 
is 9tQ[r], this entails that 



dtQ[T]=rx Q[t]. 

On the other hand, by definition of -tr[c] and iJrM, we have 

a ~ Lr[(j] + Rt[(j] X Q[t](x — xo), 
so that, differentiating in time and using (|52p. we obtain 

dt<T = dtLr[<T] + {dtRrW]) X Q[t]{x - Xo) + i?rM X (r X Q[t]{x - Xq)). 
Operating r^^^ on the right for both sides of the previous equality, we get 

i9tcr(i,Tj~^(x)) =dtLr[a] + dtRr[(^] x {x - xb) + R^'^] x (^r x (x-xb)). 
We now plug (|5T|) and (|53)) into the right hand side of ((50)) to obtain 

/ psodtT ■ dtcr dx dt ^ I (hit) + l2{t))dt, 

J[0,T]xSo ■J[0,T] 

with 

hit) := / psoiTr\x))(eit)+rit) X {x-XB{t))) ■ (OtLMit) + dtRr^W) x (x - xsit)) 

hit) := / pso 
Jsit) 

We use the identity (|22|) with 

ih, r-i, i2, ra) - (^(i), r{t), dtLr[a]{t), OtRAcrKt)), 

to get 

hit) = me{t) ■ dtLMit) + J[T{t)]r{t) ■ dtRr[a]{t). 

Finally we observe that 

hit) - / pso(r,-Hx))e{t) ■ (RAa](t) x [r{t) x {x - XB{t))])dx 
Js{t) ^ ' 



(52) 



(53) 



(54) 



dx. 



{T^^{x))(^(t)^'r{t) X (x-XB(t))) • [i?rM(<) x (r{t) x {x - xs^t)) 



dx. 



(55) 



m- (^Rr[a]{t) X (r{t) X pso{t-\x)){x - xeit)) dx])^ 



according to (211) ■ Combining (IMl), dSS]) and yields (11^1) . 



(56) 
□ 



Proof of Theorem [II 

1. Let {u,XB,r) is a classical solution of the PDEs formulation on [0,T] and let (r, 77) their respective 
flows given by formulas (|20p and (1^^ . Equation ([T]) reads 



+ (Vp) o = 0, for a; e Jo- 
Let us now consider (cr, /i) in T(^^ jj^C. Using (|57l) and performing a change of variables, we get 



(57) 



C/,;[/i] o rjdx ~ — I Vp • C/r;[M] <^2;. 
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Now using (1551) , ([55]) and Green's formula we deduce 

f d'^V f 
~ / "HIT ■ ^'fl^l ° ^'^^ = / vUj\[i\-ndV. 
JT{t) at JdS{t) 

Now using and then ©-([I]), we obtain 

~ -^-Uy^Mo-qdx / p{Lr[(j]+ Rrlcr] -X. {x - XB{t))) ■ ndT{x) 

Jj^(t) JdS{t) 

= mxUt) ■ LM + {J[Tmr{t)y ■ RM- 
Then we integrate by parts over [0,T] and conclude with Lemma [5] that (t, 77) is a geodesic on C. 

2. Conversely, let S £ be a geodesic. Using again Lemma [U it means that for any (cr, ^) G T(^ ,j)£, 
one has 



mi{t) ■ Lr[a]{t) + J[T{t)]r{t) ■ Rr[a]{t) + / dtr]{t,x) ■ ^t{Ur[^Ji]{t,r|{t,x))) dx) dt ^ Q. (58) 
[o,T] ^ Jto ' 

We first use ((58|) with ct = 0. We consider w : [0, T] x J^t — > satisfying div [w] = in J't, w ■ n — ^ 
on dTt and u'(0, •) = and w(r, •) = 0. Then (cr, /i) G T(,-^)£ for cr = and fi{t,x) :— w{t,r]{t,x)). 
Consequently, one has 

0— / dtri(t,x) ■ dt{w{t,ri{t,x))) dxdt 

df'r]{t, x) ■ w{t, r]{t, x)) dx dt 

[0,T] JFo 

dirj{t, ri^^{t, x)) ■ w{t, x) dx dt. 

It follows that d^rjit, rj^^it, •)) is a gradient field in J"t, that we denote — Vp. Now going back to ([58]) we 
get that for general (cr, /i) G T(^^)>C, one has 



m£{t) ■ dtLr[<j]{t) + J[T{t)]r{t) ■ dt{Rr[a])it) + / Vp(t, x) • C/„[Ai](t, a;) da; dt = 0. 

[0,T] ^ Jj^t ' 

Using Green's formula for the last integral and ([M)) -([55 |) . we deduce that for any cr, 

' ml{t) ■ dtLrW^) + ■ dt{Rr\o\)it) 

[0,T] \ 



+ Lr[(r]- p{t, x)n{x) dV + Rr[cr] ■ p{t, x){x - XB{t)) x n{x) dV \ dt ^ 0. 

JdTt JdTt J 

Now we integrate by parts in time the first two terms. Since this is valid for any cr, hence for any Lr[cr] 
and Rt[(t\, we infer ([3])-(|4]). Equations ([S])-® then follow from the very definition of C. □ 
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